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A quasi-spherical approximation scheme, intended to ap- 
ply to coalescing black holes, allows the waveforms of grav- 
itational radiation to be computed by integrating ordinary 
differential equations. 

04.30.-W, 04.25.-g, 04.70.Bw, 04.20.Ha 

The coalescence of binary black holes is expected to 
be one of the main astrophysical sources for upcoming 
gravitational- wave detectors. The initial phase of inspi- 
ral and the final phase of ringdown are understood in 
terms of post- Newtonian and close-limit approximations 
respectively, but the coalescence is only qualitatively un- 
derstood and generally thought to be tractable only by 
numerical methods Considerable problems have 

been encountered and currently there are no reliable pre- 
dictions of waveforms. 

This article presents an approximation scheme to com- 
pute the gravitational waveforms for space-times close to 
spherical symmetry. This is intended to apply to binary 
black holes once they have coalesced, i.e. when a marginal 
surface encloses both sources. The quasi-spherical ap- 
proximation will be best when the angular momentum is 
small, but note that even the maximally rotating Kerr 
black hole is 70% spherically symmetric according to the 
ratio of the areal and equatorial radii of the horizon. 
Thus rough estimates may still be possible even for ap- 
preciable angular momentum. 

The basic idea of a quasi-spherical approximation is 
to make a 2-1-2 decomposition of the space-time and lin- 
earize only those parts of the extrinsic curvature which 
vanish in spherical symmetry, cf. Bishop et al. Thus 
when the linearized fields vanish, spherical symmetry is 
recovered in full. This can be a highly dynamical situ- 
ation; there will be no assumption of quasi-stationarity. 
Likewise, there will be no assumption of an exactly spher- 
ical background. Unlike previous work on null-temporal 
formulations a dual- null formulation is adopted 

here, i.e. a decomposition of the space-time by two inter- 
secting foliations of null hypersurfaces. This is adapted to 
the radiation problem in that the imposition of no ingo- 
ing radiation and the extraction of the outgoing radiation 
are immediate. It also allows a remarkable simplification 
from partial to ordinary differential equations. 

A general Hamiltonian theory of dual-null dynamics 
1^ has been applied to Einstein gravity Q and is sum- 
marized as follows. Denoting the space-time metric by g 



and labelling the null hypersurfaces by x , the normal 
1-forms = —dx^ therefore satisfy 

g-i(n±,n±)=0. (1) 

The relative normalization of the null normals may be 
encoded in a function / defined by 

= -g-\n+,n-). (2) 

Then the induced metric on the transverse surfaces, the 
spatial surfaces of intersection, is found to be 

g + 2e^^ n+ (g) n~ (3) 

where (g) denotes the symmetric tensor product. The dy- 
namics is described by Lie transport along two commut- 
ing evolution vectors u±: 

K,«-]=0. (4) 

Specifically, the evolution derivatives, to be discretized 
in a numerical code, are 

A± = (5) 

where _L indicates projection by h and L denotes the Lie 
derivative. There are two shift vectors 

s± = ±u±. (6) 

In a coordinate basis u_; e^) such that u± = d/dx^, 
where = d/dx'^ is a basis for the transverse surfaces, 
the metric takes the form 

g = hij{dx' + s\dx'^ + sldx^) ® 

{dx^ + s-^(ix+ + s^_dx^ ) - 2e^^ dx+ (g) dx~ . (7) 

Then {h, /, s±) are configuration fields and the indepen- 
dent momentum fields are found to be linear combina- 
tions of 

e± = *L±*1 (8) 

a±^±L±h-9±h (9) 

i/± = L±f (10) 

u;^^efhi[l^,l+]) (11) 

where * is the Hodge operator of h and L± is shorthand 
for the Lie derivative along the null normal vectors 

l±^u±-s±=e-fg-\72^). (12) 
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Then the functions 6± are the expansions, the traceless 
biHnear forms a± are the shears, the 1-form u is the 
twist, measuring the lack of integrabihty of the normal 
space, and the functions i'± are the inaffinities, measur- 
ing the failure of the null normals to be affine. The fields 
(9± , CT-t , , a;) encode the extrinsic curvature of the dual- 
null foliation. These extrinsic fields are unique up to du- 
ality ± I— > ^ and diffeomorphisms which relabel the null 
hypersurfaces, i.e. dx^ i-^ e^^dx^ for functions \±{x^). 

It is also useful to decompose h into a conformal factor 
r and a conformal metric k by 

h = r^k (13) 

such that 

A±*1 = (14) 

where * is the Hodge operator of fc, satisfying >i=l = ir^. 
Denoting the covariant derivative of h by Z?, the Ricci 
scalar of h is found to be 

R ^ 2r-'^ {I ~ D"^ \nr) (15) 

by using the coordinate freedom on a given surface to fix 
k as the metric of a unit sphere. 

The dual-null Hamilton equations and integrabihty 
conditions for vacuum Einstein gravity have been given 
previously in a slightly different notation, so will not 
be repeated here. They are linear combinations of the 
vacuum Einstein equation and a first integral of the con- 
tracted Bianchi identity. This is the vacuum Einstein 
system in first-order dual-null form. The vacuum case 
suffices for the application, outside the black holes. 

In spherical symmetry, (s±, cr-t, Z?) vanish, while 
(h, f,d±,i/±,A±) are generally non-zero, e.g. |]lO| , p]| . 
The quasi-spherical approximation will therefore consist 
of linearizing in {s±,a±,u}, D). In practice, one truncates 
the equations by setting to zero any second-order terms 
in {s±,a±,uj, D). This greatly simplifies the equations. 



leaving the momentum definitions as 

A±r = ^re± (16) 

A±/ = i/± (17) 

A±k^r-'^(j± (18) 

A+s^ - A^s+ ^2e-fh-\uj) (19) 

and the remaining equations as 

A±0± = -v±e± - \el (20) 

A±e^ = - e-i'r-^ (21) 

A±J/^ = - e-fr^^ (22) 

A±a^ = ^{e±a^ - 9^a±) (23) 

A±w = -e±Lj ± \{Dv± - De± - e±Df). (24) 



This follows immediately from the full equations |^], the 
above expression for R and the fact that 



A± = ±i± (25) 

in this truncation. One may take quasi-spherical coor- 
dinates X* = (i?, (fi) on the transverse surfaces such that 
i 1 = sin dddAdip, the standard area form of a unit sphere. 
Then r is the quasi-spherical radius. 

The shear equations, composed into a second-order 
equation for fc, become 

□ fc = (26) 

where □ is the quasi-spherical wave operator: 

n(l> = ~2ef (A(+A_)0 + 2r-iA(+rA_)<?!)) . (27) 

Thus the conformal metric k satisfies the quasi-spherical 
wave equation. Then k may be interpreted as encoding 
the gravitational radiation. In particular, fixing to 
be the outgoing direction, the Bondi news at future null 
infinity 3+ is essentially r~^a- ||l^, as described explic- 
itly below. Likewise, the no-ingoing-radiation condition 
is just r~^a+ = at past null infinity 3^. That k gen- 
erally encodes the free gravitational data was suggested 
by d'Inverno & Stachel |13[| and has been rediscovered by 
various authors, e.g. 

The dual-null initial-data formulation is based on a 
spatial surface S and the null hypersurfaces S± locally 
generated from S in the u± directions. The struc- 
ture of the field equations shows that one may specify 
{h, f ,0±,uj) on E, ((T-(_,z^+) on (s_,cr_,j/_) on E_ 
and s+ in C/, a region to the future of E±. In partic- 
ular, the initial data is freely specifiable. There are no 
constraints as in the 3-1-1 formulation; these have been 
converted into evolution equations along I]± , which even 
in the general (non-quasi-spherical) case can be solved in 
closed form [H. 




z 

FIG. 1. The dual-null integration scheme. Initial data is 
prescribed on a spatial surface E and the null hypersurfaces 
E+ and E_ generated from it. 

A numerical integration scheme runs as follows, as de- 
picted in Fig.|l|. First integrate the A+ equations from 
E to obtain the full data on E+. Then integrate the 
A_ equations one step along each ingoing null hyper- 
surface, generating a new null hypersurface Then 
repeat: integrate the A+ equations along S'^ to obtain 
the full data on E^, and so on. In practice, some inter- 
polation between the two integrations is useful. There 
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are many ways to perform the integrations in a differ- 
ent order, allowing flexibility which can be used, for in- 
stance, to avoid singularities. Any such scheme gives two 
estimates of {r, k, f,9±,uj) at each point, since some of 
the equations play the role of integrability conditions. 
Thus one could ignore such equations to obtain a free 
rather than constrained integration scheme. This allows 
numerous internal checks on the accuracy of the numeri- 
cal code, analogous to those of 3-1-1 integration schemes, 
e.g. Choptuik |15| . 

The equations for A±(r, f,9±,v±), the quasi-spherical 
equations, decouple from the remaining equations. Thus 
there is a quasi-spherical background which may be found 
by integrating the quasi-spherical initial data. Since this 
background is independent of the linearized part, one 
may economize when computing different evolutions on 
the same background. It should be stressed that the 
quasi-spherical background is neither fixed in advance 
nor necessarily spherical, e.g. Dr in general. 

To compute the outgoing radiation, one now needs 
only to integrate the equations for (Aj-fc, Aj-CTzp), i.e. 
the quasi-spherical wave equation for k. It is remark- 
able that this entire integration scheme involves only or- 
dinary differential equations. The equations for A-tcj 
are partial differential equations, containing transverse 
D derivatives, but the other equations decouple from the 
equations for (s±,ti;), which therefore need not be solved 
for the radiation problem. In short, most of the complex- 
ity of the system has been isolated and sidestepped. 

Moreover, one may use a conformal transformation to 
obtain a scheme which is more accurate at large dis- 
tances. Using the conformal factor 



the rescaled expansions and shears 



§± ^ r0± 



r ^a± 



(28) 



(29) 
(30) 



are finite and generally non-zero at 5^ for an asymptot- 
ically fiat space-time. Rewriting the relevant equations 
yields the quasi-spherical equations 



A±z?± 



and the linearized equations 

A±fc = n<;± 



(31) 
(32) 

(33) 
(34) 

(35) 



(36) 
(37) 



One may take S+ to be either part of 3 , as depicted 
in Fig.H, or at sufficiently large distance for numerical 



purposes. Here, large distance means small ft. For the 
quasi-spherical approximation to be valid at large dis- 
tance, one may fix {f,'&±,k) = (0, ±-\/2, e) on E, where 
e = d-d dd + sin^ "ddip (S) dip is the standard metric of a 
unit sphere, and = on S+. The remaining coordi- 
nate data is given by on the ingoing null hypersurface 
which is left free so that one may adapt the foliation 
of I]_ to the surfaces which are most spherical. 

The no-ingoing-radiation condition is "j-j. = at 
leaving the gravitational initial data as on S_. The 
outgoing radiation is found by computing at 3"*", 
which is essentially the Bondi news. More precisely, the 
Bondi energy fiux at 3^ would be [p2l 



5/^^ 1 1 



^±1 



647r 



(38) 



where \ \a\ 



k k GacCfbd and such second-order terms 



are no longer being ignored. That is, the energy supply 
would be 



A±E = f iip± 



(39) 



where E is the Bondi energy. In summary, the outgoing 
waveforms and their energy may be computed by inte- 
grating 9 first-order ordinary differential equations and 
their duals, or a subset in the case of free evolution. For 
numerical purposes, this is a dramatic simplification. Nu- 
merical implementation of this scheme is in progress [ p^ . 

Before concluding, it should be noted that the domain 
of validity of the quasi-spherical approximation is not 
known in a precise sense. The guarantee is simply that 
spherically symmetric Einstein gravity is recovered in full 
when the linearized fields vanish. For the usual pertur- 
bative approximations, one may check successive orders 
of approximation to compare accuracy, but for the quasi- 
spherical approximation, the corresponding second-order 
approximation would be full Einstein gravity. If one 
wishes to know whether a given space-time is sufficiently 
spherical, the rough answer is that there should be a 
2-1-2 decomposition such that the fields to be linearized 
are small compared to the remaining part. This depends 
on the choice of transverse surfaces, so that there will 
be some art to choosing the 2-1-2 foliation for optimal 
accuracy. For the Kerr black hole, the quasi-spherical 
null coordinates of Pretorius & Israel may be useful. 
For a coalesced black hole, one might base the foliation 
on the marginal surfaces which locally define it, i.e. use 
the coordinate freedom in on I]_ so that the folia- 
tion contains a marginal surface. There are some general 
laws of black-hole dynamics in terms of marginal surfaces 
and the trapping horizons they generate [lqjl9[| , includ- 
ing that outer trapping horizons are achronal and there- 
fore cannot causally influence 3+. Thus S- need not 
extend inside the trapped region in order for the domain 
U of integration to reach all of 3+ . 

To conclude, the intended scenario is an asymptotically 
flat space-time containing coalescing black holes, with an 
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ingoing null hypersurface S_ chosen to intersect the coa- 
lesced black hole, i.e. the region of future trapped surfaces 
enclosing the original black holes. The initial data on I]_ 
may be determined by extracting the relevant data from 
a conventional 3+1 numerical computation from an ini- 
tial spatial hypersurface Soj smoothed off to the past of 
So- The smoothing may be expected not to affect the 
results significantly if S_ n Sq is sufficiently outside the 
black holes; there will be some spurious radiation at 9"'" 
at early times, but not at the relevant late times, as this 
would involve backscatter of backscatter. 




FIG. 2. Application to coalescing black holes. Shading in- 
dicates regions containing trapped surfaces, with the outer- 
most trapped region being that of the coalesced black hole 
and the inner trapped region that of one of the original black 
holes. Initial data on the ingoing null hypersurface E_ may 
be extracted from a conventional 3-1-1 code based on an initial 
spatial hypersurface Eo with outer boundary B. 

The advantage of this procedure is that it avoids the 
outer-boundary problems which plague the conventional 
codes. As depicted in Fig.||, the outer boundary B can- 
not causally influence I]_ if it intersects Eg inside B. 
Thus the scheme requires only clean data from the 3-1-1 
computation, uncontaminated by outer-boundary prob- 
lems. A code implementing the scheme may be regarded 
as a black box which, taking input from any other code 
from which the required data on 5]_ can be extracted, 
computes approximate waveforms for the gravitational 
radiation. 

This suggests a quite general proposal to compute out- 
going gravitational radiation from a 3-1-1 computation by 
a conformal dual-null code which extracts data on an in- 
going null hypersurface intersecting the initial spatial hy- 
persurface inside its outer boundary. One might expect 
this to be simpler than the usual matching on a tempo- 
ral hypersurface 1^-0] since the outer boundary is avoided 
and the problem is merely of extraction rather than dy- 
namic matching. At present, there seems to be neither 
a general conformal dual-null code nor work on data ex- 
traction on a null hypersurface, though the null-temporal 
formulation can presumably be adapted 0|. These are 
tractable projects which would allow accurate compu- 
tation of gravitational waveforms from coalescing black 
holes. 
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